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Abstract 

We study the matrix factorization problem associated with an 50(2) spinning top by 
using the algebro-geometric approach. We derive the explicit expressions in terms of 
Riemann theta functions and discus some related problems including a non-compact 
extension and the case when the Lax matrix contains higher-order powers of the spec- 
tral parameter. 
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1 Introduction 



Analytic factorization for matrix valued functions, or the matrix Riemann-Hilbert prob- 
lem, appears in problems of exactly integrable dynamical systems as well as in prob- 
lems of elasticity and diffraction. 

In the context of exactly integrable models, the factorization problem for the R- 
matrix associated to an affine Lie algebra is the same as the matrix Riemann-Hilbert 
problem [JTp . In this case exists a general formula for the solution in terms of Riemann 
theta functions [|TJ], derived from algebro-geometric methods of finite-band integration 
theory [g]. This construction is based on the fact that the time evolution linearizes 
on the Jacobian of the spectral curve, and consequently one can express the solution 
in terms of the Baker- Akheizer functions, which are eigen-vectors of the Lax matrix. 
Various integrable dynamical systems can be formulated via such R- matrix [Jl[], and 
most typical are the spinning tops. 

In this paper we derive the explicit expressions for the case of an 5*0(2) spinning 
top, which is a simple but a non-trivial example of the two-by-two matrix factorization 
problem. This is done in order to illustrate the general theta-function formula from [p, 
and also in order to be able to compare to alternative methods of matrix factorization. 
The neccesity for alternative methods of factorization comes from the fact that the 
theta-function formula becomes difficult to use in the case when the spectral curve is 
non-hyperelliptic, which is mostly the case. Also there are cases when only a non- 
canonical factorization is possible, for which the theta-function formula is not valid. 
We illustrate this on the example of a non-compact extension of the 5*0(2) spinning 
top. 



2 Lax formulation and matrix Rieman-Hilbert prob- 
lem 

Consider a rotator with a unit radius and angle in a potential V(<fr) = —a 2 cos 20, 
where a is a real constant. The Hamiltonian is given by 

H = h 2 -a 2 cos 20 , (1) 

where I = <p = d(f)/dt. As a one-dimensional Hamiltonian system, it is integrable, and 
it is the n = 2 case of an n-dimensional integrable spinning top system given by the 
following n-by-n Lax matrix ^ 

L = aA + 1 + sA" 1 = aA + 1 + R^fR^X' 1 , (2) 

where A is a complex number called the spectral parameter, a and f are fixed symmet- 
ric matrices, 1 is antisymmetric and traceless matrix, while is an 50 (n) matrix. 
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Formula (g) gives a mapping from the phase space of the n-dimensional spinning top 
T*SO(n) = (Rfa 1) to a Hamiltonian orbit in the affine gl(n) algebra g defined as 

9 = Y1 g 'i XJ ' 9j £ gK n ) ■ 
jez 

Without loss of generality one can take that the a and f matrices are diagonal, and for 
simplicity we take a = f and tra. = 0. The last condition implies tr L = 0, and this 
means that we are restricting g to the affine sl(n) sub-algebra. 

The Hamiltonian ([]]) is the n = 2 case of the SO(ri) Hamiltonian 

H = ~Res tr(\- l L 2 ) = -hrl 2 - tr{aR^HR^) . (3) 

The corresponding Lax equations follow from the R-matrix formulation [[]]], and they 
are given by 

L=[L,M±] , M ± =~(R±l)d<p(L) = ±P ± L , (4) 

where (p(L) = \tr(L 2 ) is an invariant function on g and R — P + — P_ is the R-matrix. 
P± are the projectors onto g± subalgebras, which are defined as 

9+ = ^9j^ , fj- = ^2gj\ 3 , (5) 

j>0 j<0 

so that g = g + © g_, and therefore 

M+ = P+L = 1 + aA , M_ = -P_L = -sA" 1 . (6) 
The equations (ffl) are solved by 

L(\,t)=g ± (X,t)L(\)gg\\,t) (7) 

where L(X) = L(\,t = 0) and g±(\,t) are solutions of the following matrix factor- 
ization problem 

e tL ^ = g- l (\,t)g4\,t) , (8) 

such that g + is analytic for A ^ oo and g_ is analytic for A ^ 0. This is a Riemann- 
Hilbert problem for e tL ^ and a curve which encircles the A = point in the complex 
plain. 

Solution of the factorization problem (|8|) exists for t that satisfies the Gohberg- 
Feldman bound, which is satisfied for sufficiently small times O. Otherwise one has 
a non-canonical factorization 

e tL ^ = g- 1 (X,t)D(X)g4X,t) , (9) 

where D(X) = diag(X kl , A fcn ). 
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3 Algebro-geometric factorization 



One way to solve the problem (^]) is to use the connection with Riemann surfaces. The 
spectral curve is given by 

det(L(A) - y) = , (10) 

and it is time-independent. Its compact model defines the corresponding Riemann 
surface T. Consider the line bundle E L (p) e C n associated to the eigen- vectors of L 

L(p)v(p) = fj,(p)v(p) , per. (11) 

Its time-evolution is linear on JacT, and it is given by E t = E ® F t , where F t is the 
line bundle on T determined by the transition function exptfi(p) with respect to the 
covering {U+, LL}, where U± = {p \ X ±1 (p) ^ oo}. In order to solve (||), one needs 
to consider the holomorphic sections of the dual bundle E L , since the bundle El does 
not have holomorphic sections. Denote these sections as ip(p), then they satisfy 

L(p)ip(p) = fi(p)ip(p) , (12) 

so that ip±(p, t) = g±(X(p),t)i/j(p) are eigenvectors of L(A, t) which are regular in U±. 
Hence 

g ± (\,t) = ^ ± (X,t)^l 1 (X,t = 0) (13) 

where \1/±(A, t)jk = ^(Pk, t). Hence if one knows the sections tp±, one can solve the 
factorization problem @ as well as the dynamics. These sections can be represented 
via meromorphic functions (p± on U± as 

<P±(p,t) = , (14) 

where s±(p) are sections of E* L . These functions are called Baker- Akheizer (BA) 
functions and satisfy 

(<P±) > ~D* (15) 

where D* is a divisor of degree g + n — 1, where g is the genus of V. 
Given ( |T3[ ) and (JT3J) and the following properties of ip± 

d ^L = -M ± (X,t)^,i) (16) 

and 

^ + (p,t) = e-^V-(P,*) , d7) 
one can construct the BA functions in terms of the theta functions []3], |], [TJ]. Let = 
Y^j=i Pj t> e me divisor of poles of A, let po be a fixed point in T and let D be a divisor 
of degree g such that the line bundle El is associated to the divisor 

D* = P oo + D-p . (18) 
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Let Uj be a set of normalized holomorphic Abelian differentials on V 



u k = 2ni5 jk , uj k = B jk , (19) 

a,j Jbj 

where aj and bj, j = 1, g is a basis of homology cycles and £> is the period matrix. 
Note that the Riemann theta function is defined on C 9 by quasi-periodicity conditions 

9{z + 2mn) = 9{z) , 9{z + Bn) = e -^ Bn ' n) - {z ' n) e(z) , n e Z 9 (20) 
and can be written as 

0(s) = ^ exp{^(Bn,n) + (z s n)} , (21) 
nezs 

where (X,Y) = Y^j=\XjYj. The theta function can be considered as a section on 
JacT = C ff /(27rm + Bm). 

Let A(p) be the Abel transform 

A(p)= f uj , (22) 

and let d£l be an Abelian differential of the second kind, normalized by 

dVL = , (23) 

which is regular in C7_ such that dfi — d/i is regular in U + . Then 



^ fl, t) - 7 f^ fa ) ^(P)-C-^)^(P)-Q)) 



(24) 



where = f p dQ, and the constants c k and the function 7 3 -(i) are determined from 
the following requirements 

(a) has no discontinuity across the cuts 

(b) Lp 3 _ is subordinate to the divisor p — Pj — D 

(c) itesp.(A-yi) = dj 

where dj are time-independent constants. 

Condition (a) is satisfied if c + c = c\ + c, and Vj ; = J b dfl. Condition (b) is 
satisfied if the constants c k are chosen such that the divisor of zeros of 9(A(p) — Co) 
is po + D, the divisor of zeros of 9(A(p) — C\) is D + D and the divisor of zeros of 
9(A(p) — Cj) is Pj + D, for some divisor D. Condition (c) serves to determine the 
functions jj(t) and it follows from dl6| ) when A — > oo. The following identity will be 
useful for our purposes 

» n 

^ = / dtt= - ^Res Pk {pujj h ) . (25) 
k=i 
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4 The SO(2) Case 

In the case of the rotator, a = f = diag(a, —a), 



1 ( ° I ) • (26) 



and 



so that 



-I 



cos 20 sin 20 
sin 20 — cos 2<? 



(27) 



t — ( a ^ + a ^ 1 cos 20 I + aX 1 sin 20 \ 

\ -/ + aA" 1 sin20 -aA - aA" 1 cos 20 J " ( 8) 

The spectral curve ( JT0| ) is then given by 

= a 2 A 2 - 2E + a 2 A~ 2 , (29) 



where E = H is the energy. The equation (29) can be transformed into elliptic form 



w 2 = A 4 -^|a 2 + 1 = (A 2 -/ 2 )(A 2 -/ 2 ) , (30) 

gt 

so that the spectral surface is a torus. We take the a± cycle to encircle the [— 1 + , — /_] 
cut, while the b\ cycle goes from the [—1+, — /_] cut to the [/_, l + ] cut, were it passes 
through to the second sheet and it goes back to the [— /+, — 1-] cut. The normalized 
holomorphic Abelian differential is given by 

2ni dX f dX ,„.,. 

w = -: , A= — , (31) 

^ w J ai w 



so that the period matrix is given by 

r b 

B — I lu = 2m— = 2mr . (32) 

Jbi -A 

The Abelian differential of the second kind dtt must satisfy for A — ► 

a 



dn( P ) = dfi + o(i)dx= =f— + o(i) 

A z 



rfA , (33) 



and it is normalized by 

' dtt = . (34) 



Hence the second-order poles are at p = + and p = CT and one can construct dVl in 
terms of the prime forms 0]. Alternatively, from ( |33| ) it follows that the corresponding 
Abelian integral fl(p) behaves as 

n(p) = ±~ + 0(l) , (35) 



for p — > C^, and hence it is a meromorphic function with simple poles at (J 1 * 1 , and with 
residues ±a. A meromorphic function having only simple poles can be specified by 
the position of the poles and with the values of the residues, and it can be constructed 
in terms of the theta functions This implies that 



fl(p) = aa 



9[(A(p) - A(0+)) 9[(A(p) - A(0-)) 
e^Aip) - A(0+)) 9 x {A{p) - A(0r)) 



(36) 



where a = 2ni/ A and 9\{z) is the odd Jacobi theta function ( p0| ) and 9'{z) - — — 
Also note that the identity (E51) gives 



dz 



2ni 

V = 2a — = 2aa . (37) 

We take D = p\ ^ e, ; , where are the end-points of the brunch-cuts, and P^ = 

P x + P 2 — oo + + oo~ so that 



and hence 



By using 



where 9 x (z) 



c = m + -B 

C X = A{ Pl )+7Ti + ^B 

Cj = A{P-)+m + ^B 



c = a + cj - cq = A(p x ) + A(Pj) + ni + -B 



(z + m + -B) 



-B/S-z/2, 



-i9[\, \}{z) is the odd Jacobi theta function, we obtain 



9 1 (-A(p) + A( Pl ) + A{P 3 ) + Vt^-Ajp)) 
9 X {-A(p) + A(p x ))9 x (—A(p) + A(P,)) 



(38) 



(39) 
(40) 

(41) 



and Cl(p) = fi(p) - V/2. 

The constants jj(t) can be determined from the condition (c), and for this one 
needs the expansion of ( |4"T| ) in ( = A -1 for A — > oo. The following formulas will be 
useful 

4- 

cp 



A(p) 



pi) 



LU+ LU = A(P j 



and 



« i c + o(c 2 ) = ±«c + o(C 2 ) . 



One also has for p — > P~ 



n(p) = n° 1 + n)( + o(( 2 



(42) 
(43) 

(44) 
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so that when p — > P 3 - we obtain 



c 1 + 
*?(o) 



a,- 0i(Ai + yt) 



0(C) 



0i ( — A?) ^(-Aj + Ai) 2^(0), 
where A(pfe) = A^. Also when p — > P*, where p* = (—w, A), we have 

,oo 9i(A! + A* - A* + Vt)6A-A*:) 



(45) 



6 1 (-A* j +A 1 )e 1 {-A*+A j ) 



(46) 



From (03]) it follows that 



(47) 



9 1 (-A ] + A 1 )a J 9'M ' 
As far as the solutions of the equations of motion are concerned, these can be 



obtained by inserting the A — > oo expansions ( Hq ) and ( ftq ) of the BA functions into 
(fl6|). Namely, for A — > oo we have 



so that (16) gives 



^ = A^i = Xipi + i){ + 0(1/ X) 



(48) 



(49) 



By using (g9|) together with (|45j) and @ we get ^ = 0) T , = (°> rf 2 ) T and 



4- rf 
± a — 

dt 



_e 1 {A 1 + vt) 



d_ 

dt 



• ^ g 1 (A J -^ + A 1 + W)g 1 (-A*) 
7,uc 3 eMi ~ ^MAj - A*) 



±a cos 20 
adj sin 20 



(50) 
(51) 



From (50) we obtain 



sm 



-I (ni /a-l)-^^log^) 



(52) 



where Q t = Q\ = -fif and z = Vt + Ai. Due to (|23J), V^/2a = a so that (g2|) can be 
written as 



sin 2 = const. + a 2 p{2aat + A\) , 
where p is the Weierstrass function. Since 



(53) 



(54) 
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where u>i 2 are the periods of the Weierstrass function, we get 



sin 2 <j> = p ( 2at H 1 + const. . (55) 



a 



The expression (55) can be shown to coincide with the solution of the equations of 



motion obtained from the energy expression 



2E = 2 -2a 2 cos 20 . (56) 



It follows from (|56|) that 



ut = ± I , (57) 

Jo y 1 — k 2 sin 2 

where Q = y/2(E + a 2 ), /c 2 = 2a 2 (£ + a 2 )- 1 and E > -a 2 . The elliptic integral (|57|) 
implies 

sin = sn(u;t, /c) . (58) 



However, in order to compare to (]55|), we rewrite the integral (|57|) as 



siu- 



dx 

2at = ± I ^^^=^^^= . (59) 

y4a;(l — x){k 2 — x) 



o 



Since the Weierstrass function can be defined as 



/i3 dv 

Jp{u) V% - 92V ~ 93 

from ( |59"| ) it follows that 

sin 2 <^ = <p{2at + -u ) + const. , imw ^ , (61) 

where 

r°° dx 

u = / . . (62) 

Jo ^Ax{l - x){k- 2 - x) 

Clearly the expression d6ll ) coincides with (|55|). 

Now define the matrices $±(A, t)jfc = <p±(pk, t), where 

^ + (p,t) = e-^V-(P,*) > 

then 

5± (A,t) = $ ± (A,t)$± 1 (A,t = 0) . (63) 
This formula follows from (jl3|), due to relation (jl4|). 



8 



The relation ( |5"5l ) can be satisfied for all times if the choice of p\ is such that 
Im{Ai/a) 7^ 0, in accordance with (|6T|). This follows from p(u) = u~ 2 + 0(u 2 ) 
so that the r.h.s. of ( |53| ) diverges for 



2at + — = . (64) 

a 

Therefore if Im(Ai/a) ^ 0, one obtains a physical solution. 

If pi is chosen such that Im(Ai/a) = 0, then the factorization (f|) will not be valid 
for t close to A\/(2aot). In this case one would have a non-canonical factorization 
given by @. We will show in the next section that this case corresponds to a dynamical 
system based on a non-compact extension of the 5*0(2) group. 



5 Related problems 

Let us replace </> with iq in ([I]), and let I = q. In this way we obtain a new dynamical 
system describing particle on a line moving in a potential V(q) = —a 2 cosh2g. The 
solution of the equations of motion can be obtained from the energy 

2E = q 2 - 2a 2 cosh 2q , 

and it is given by 

r inh2 q dr 
2at = ± / , (65) 

Je k a/4x(1 + x)(k- 2 + x) 

where e k = for k 2 = 2a 2 /{E + a 2 ) > 0, while e k = -&T 2 for k 2 < 0. Note that 
now there is no restriction E + a 2 > as in the compact case, since the potential is not 
bounded from bellow. Then by using the same relations as in the compact case we get 

sinh 2 q = p(2at + u ) + const. , Imu = , (66) 

since 

f°° dr 
u = / = . (67) 

Je k V /4X ( 1 +X){k~ 2 + X) 

This is an example of a dynamical system where q diverges for finite times, and 
hence one can expect that the corresponding factorization will be non-canonical. The 
replacement <p — > iq corresponds to replacing the 50(2) group with the non-compact 
subgroup of 50(2, C). The corresponding Lax matrix can be obtained from (0) by 
taking f = —a and <\> = iq, so that 



t_( aA — aA 1 cosh 2q I — ia\ 1 sinh 2q 
\ — I — iaX^ 1 sinh 2q — aX + aA -1 cosh 2q 



(68) 
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and I = iq. The spectral curve is ( [Z9| ) with E — > —E, since 

-I 2 - 2a 2 cos 20 = q 2 - 2a 2 cosh 2g = 2E . 
Consequently, the derivation of the BA functions is almost the same as in the section 



4, and the only difference is in the equation (pOj), where a — > —a due to the fact that 
now f = —a. This gives 



sinh 2 q = p(2at + Ai/a) + const. . (69) 



The physical solutions are obtained for Im Ax/a = 0, in agreement with dog). Hence 
the canonical factorization of e tL will not be possible for 2at — > Ax/a. 

Note that for n > 3 one obtains the spectral curves which are non-hyperelliptic, 
and in that case finding the BA functions gets difficult, because not much is known 
about the non-hyperelliptic curves, except in some special cases, like Kowalewski top 
0] ( see [0 for other examples). For example, for n = 3 case the Lax matrix (Q) gives 
a cubic spectral curve 

/i 3 + L (2) /i - det L = , (70) 

which is difficult to analyze. It would be interesting to see whether some alternative 
method of solving the corresponding factorization problem can be used, for example 
see [§]. Preliminary study indicates that the techniques developed in [j^] can be used 
to solve the n = 2 case [[7p, so that one can expect that the n = 3 case could be also 
solved. 

Note that the Lax matrix (0) can be generalized to 

-fc 

L = a\ m + L i xj > m,keN , (71) 

j=m-l 

where a is a constant matrix while Lj are the dynamical matrices [JTp . The dynamics is 
given by (Q) and the corresponding factorization problem is given by (Q). In the n = 2 
case (L e gl(2) or sl(2)) the BA functions can be easily calculated since in that case 
one obtains a hyper-elliptic spectral curve 

w 2 = P 2m+2k (\) , (72) 

where P2m+2fc(A) is a polynomial of order 2m + 2k. However, this system is not 
physically interesting, i.e. the corresponding Hamiltonian 

H = —-Restr(X~ 1 L 2 ) , (73) 
2 

has no physical applications. Still, the two-by-two matrices ( |7T| ) provide non-trivial 
solvable examples of the Riemann-Hilbert problem for e tL , and one can compare the 
results to alternative methods of matrix factorization. 
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